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Abstract 

In the hadronic fireball phenomenology of Gamma Ray Bursts (GRBs), it is ex- 
pected that the observed photons are accompanied by UHE neutrinos, which have 
not been observed yet. It is one of the challenges of experimental UHE neutrino 
astrophysics to look for a signal from GRBs. In this paper, the differences be- 
tween a search for a diffuse signal and an examination of a source sample given by 
e.g. BATSE will be analyzed. Since redshift information is needed to determine the 
correct energy spectrum, long duration bursts with redshifts from different estimate 
methods will be used. We will start with an overview of the current understanding 
of GRB neutrino physics and will then use this knowledge to make predictions for 
a coincidence flux and a corresponding diffuse flux. It can be shown that shape and 
normalization of the spectrum is highly dependent on the set of bursts used and 
that individual bursts can determine the total spectrum. 
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1 Introduction 



The prompt GRB photon spectrum, A^^, is usually given by a Band func- 
tion (|l|) which can be approximated by a broken power law, 



E^-< for Ey>e\. 



The spectrum is presumably produced by synchrotron radiation of electrons 
in the internal shock fronts of the jet, see e.g. There are two approaches to 
explain the break in the spectrum at a break energy of typically ~ 250 keV: 
The most common explanation is the steepening of the spectrum by a power 
of one due to cooling of electrons at high energies, see for example (0; 0) as 
a review. The break can, however, also be explained by assuming an Inverse 
Compton scattering scenario, see e.g. (j^ and references therein. Throughout 
the paper, all given energies are in the observer's frame at Earth unless de- 
clared otherwise. The spectral indices are usually scattered around average 
values of ~ —1 and ~ —2. Assuming hadronic acceleration in the jet, a 
prompt neutrino flux that is correlated to the photon spectrum results from 
photohadronic interactions in the source. The neutrino spectrum, dN^/dE^, 
can be derived assuming that the proton spectrum follows the electron spec- 
trum of the source. Since the neutrino flux in turn follows the proton spectrum 
in first order approximation, it can thus be connected to the observed syn- 
chrotron spectrum of the sources and can be described as 

(E./e^)-"" for E, < ei 

{Ejel)-^'^ for ei<E,< el (2) 

{Ejelr^'^iEjel)-^ E, > el . 

V 
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The photon spectral indices can be used to describe a,y = (3j + 1 and (3^ = 
+ 1. The second break at E,y = el in the neutrino spectrum results from the 
fact that pions lose energy at very high energy due to synchrotron radiation. 
Thus, less neutrinos result from pions at very high energies which leads to 
a steepening of the spectrum by a power of 2. A detailed derivation of the 
neutrino spectrum as presented above is given in (j^. 

The spectrum is normalized to the 7-ray fluence which is assumed to be 
proportional to the neutrino luminosity, 

x-F^= I -^dE,^ 1x1(10) -A,. (3) 
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All parameters occurring in following calculations are listed in table 1. The 
factor X is given by the product of the fraction of proton energy transfered to 
the pions, /tt, a factor 1/8 since half of the photohadronic interactions result 
in four neutrinos and a factor 1/ fe to account for the fraction of total energy 
in electrons compared to protons in the jet (0). The normalization constant 
is therefore given as 

11 F 



In the following, the normalization of a single burst will be modified to a 
quasi-diffuse normalization by multiplying A^ with the number of bursts per 
year (2/3 • 1000 long duration bursts per year) and dividing the result by 4 7r 
sr, 

3 yr 4 TT sr 



The first break energy in the spectrum, e^, is related to the break energy in 
the photon spectrum as 

{ml - ml) ■ . .-1 



It is determined through the minimal energy necessary to produce a A-resonance 
in the shock fronts of the bursts. Using the numerical values given in (0) for 
the proton mass, = 0.94 GeV, and the A mass, m^ = 1.23 GeV, leads to 

e^ = 7-10^-(l + 2)"2Jl^ GeV. (7) 

^■y,MeV 



The second break energy is connected to the synchrotron loss time. It depends 
on the neutrino flavor and for muon neutrinos, it is given as 



3 TT ee c ty 



Ar^aTesL^ {I + z) ■ me 



For electron and anti-muon neutrinos the break energy is about an order of 
magnitude lower, since these neutrinos result from the muon decay. The muon 
lifetime is about a factor of 100 higher than the pion lifetime, which reduces the 
threshold of synchrotron losses. The derivation of the second break energy can 
be found in (j7|). Here, and eg are the fraction of the burst's internal energy 
going into the magnetic field respectively into electrons. The equipartition 
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fractions have been set to eg = 0.1 and eb = 0.1. There is no good way of 
determining the equipartition fractions theoretically yet. However, afterglow 
observations indicate values of the order of 0.1 (0). The remaining parameters 
in Equ. (8) are listed in table 1 with the values as used in the following 
calculations. Inserting all numerical values gives 

< = -^/J' ^t'' yZfGeV . (9) 



The boost factor F is constrained to 100 < F < 1000, since for boost fac- 
tors less than 100, the medium would be optically thick to photons and for 
F > 1000, since protons lose most of their energy to synchrotron radiation. 
The possibility of fluctuating F using the photon break energy is given as 
demonstrated in (0), but there are several arguments for using a constant 
value: Bursts can be misaligned which would lead to a misinterpretation of 
the boost factor. Also, varying the break energy for each single burst might 
implicitly already include boost factor fluctuations. Thus, a constant boost 
factor of F = 300 is used in following calculations. 



In previous papers, e.g. (|7|), a variation of the fraction of energy going into 
pions has been discussed. Such a variation would further increase the width 
of the distribution of the neutrino spectrum normalization. varies with the 
burst luminosity, the boost factor, the photon break energy and the variability 
time as 

^ . (10) 

2.5 (^-y^MeV 



In the following, this fraction will be kept at a constant value of f^^ ~ 0.2 for 
the following reasons: 

• /tt strongly depends on the boost factor F which will be used as a constant 
as discussed before. The dependence on the other three parameters is only 
linearly and thus not as striking as a variation of F would be. 

• The main uncertainties in the current calculations result from the lack of 
knowledge of the parameters eg and as well as from uncertainties in the 
redshift relation which leads to uncertainties in L^. These three parameters 
are all important in the determination of the spectrum normalization and 
thus, a constant value is favorable. 

The variability time appears in Equ. (9) as well. Apart from the arguments 
mentioned for the normalization of the spectrum already, the effect of vary- 
ing here would not be very significant in the detection rates, because the 
dominant contribution seen in the detector comes from the energy region of 
the first break. 
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The diffuse neutrino ffux prediction derived from the model above is given in 
1^, where the authors use average parameters to determine the shape of the 
spectrum. Cosmological evolution of the sources have been considered in that 
model by taking into account the redshift evolution of GRBs. It is assumed 
that GRBs follow Star Formation Rate, since they appear to be connected to 
Supernova-Ic explosions. The final result of this work will be compared to this 
standard flux. To do this, the diffuse Waxman-Bahcall flux has to be weighted 
by a factor 2/3, since this work only deals with long duration bursts which 
make up approximately 2/3 of the total burst population. 

A different approach to predicting the neutrino flux from GRBs is to look at 
each burst individually and add up the individual spectra to make a prediction 
of the total flux from these sources. An approach of using individual spectra 
as discussed by Guetta et al. (0) will be used here to get an estimate of the 
flux that can be expected to be observed by coincident measurements of the 
AMANDA ^ telescope. In addition, a diffuse prediction will be made, using the 
mean values of the parameter distributions that are given or can be derived 
from BATSE data. 

Two different burst samples based on the determination of burst redshifts 
using the redshift estimators variability (568 bursts) and lag (292 bursts) will 
be examined in this paper. A more detailed discussion of the sample properties 
will be discussed in section 2. 



2 The burst samples 



Due to the necessity of detailed optical afterglow observations for the deter- 
mination of GRB redshifts, there are only about 61 bursts (by Oct 2005) with 
observed redshifts today. A histogram of the 61 observed redshifts is presented 
in Fig. 1. All 61 redshifts observed until today are shown as the solid line, the 
dashed line represents the sample of GRBs with redshift until June 2000, the 
BATSE era. Because of the low statistics of bursts with directly observed 
redshifts, different redshift estimators have been developed. Note that these 
methods apply only to long duration bursts. Although first measurements of a 
short duration burst redshifts succeeded recently, see e.g. 
not enough statistics for short bursts yet. Two burst samples will be used in 
the following using different estimator methods for the redshift determination. 
Both burst sets are subsamples of the BATSE catalog: 

A subsample of 568 bursts with redshifts determined using the variability of 
GRBs to estimate redshifts is given by Guetta et al. (0)- Throughout the paper, 

^ Antarctic Muon And Neutrino Detector Array 
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Fig. 1. Redshift distribution of 61 GRBs with measured redshifts (1997-Oct 2005, 
solid line). During the operation of BATSE (<June 2000), 20 bursts (dashed line) 
have been assigned redshifts. 

it will be referred to the variability sample when this set of bursts is discussed. 
Additionally, a sample of 292 bursts with redshifts determined connecting 
tem por al lag between the signal in BATSE's different energy channels is given 
by (flih . based on the work (14) and (jish . Band fits have been applied to fit 



these spectra. We will refer to this set of GRBs as the lag sample. 

The possibility of analyzing bursts with respect to their neutrino signal is 
only given for a small set of bursts because of constraints from the running 
UHE neutrino experiments. During BATSE's operation time, the AMANDA 
experiment is one of the most interesting detectors for GRB analysis: There 
are 105 bursts in the BATSE catalog that can potentially be examined using 
the AMANDA experiment now integrated into the IceCube detector as it is 
pointed out in (Il6l ). None of these bursts is given in the lag sample. There 
are however 82 bursts of the variability sample which are in the set of the 
105 bursts as well. These bursts will be discussed separately in this section. 
In the following, we will refer to these bursts as the AMANDA suhsample for 
simplicity. A work particularly dedicated to these about 100 of these bursts 
will be presented soon in the context of the analysis of a potential neutrino 



signal from these bursts ()17[ ). Not all of the 105 bursts can actually be used 
for an analysis due to the restricted availability of data on either AMANDA's 
or Batse's side. A pioneer analysis of a potential neutrino signal from monster 
burst GRB030329 in AMANDA is described in 

In this section, the parameter distributions will be discussed. Each individual 
set of bursts will be used in section 3 to make a prediction of the flux from 
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the given source sample. Additionally, average values will be determined by 
calculating the weighted mean value of each distribution as 



P = 



(11) 



with IcT errors 



a{p) = {pf - p2 



(12) 



Here, Pi is the corresponding parameter value of a certain bin i and w{pi) is 
the number of entries in that bin. n is the total number of bins. In each of the 
following parameter distributions, the average values are indicated as vertical 
lines. The horizontal lines show the la deviation from the average. The set of 
average parameters of each sample are summarized in table 2. 

It should be noted that fit parameters as the spectral indices of the photon 
spectra, and as well as the break energy in the photon spectrum, E^, are 
not identical for identical bursts in the different samples. The reason is here 
that the three parameters are correlated. Fixing the break energy for example 
at a certain energy would influence the values of the spectral indices etc. 

2.1 Lag and Variability samples 

The expected neutrino spectrum from individual bursts can be estimated from 
the parameters of the photon spectrum as it has been discussed in section 
1. The distribution of the neutrino spectra's indices follow the one of the 
photon spectrum with = [3^ + 1 and (3^ = + 1. The distributions of 
the spectral indices and (3^ are displayed in Fig. 2 and 3. In all following 
figures, the solid line represents the distribution of the variability sample, 
while the dashed line shows the lag sample distribution. The distributions of 
the first spectral index differs between the two samples: While the lag sample 
shows a Gaussian behavior, the variability sample distribution does not seem 
to give a particular pattern. The distribution seems to be more randomly 
scattered. The mean values for the first spectral index are al""^ = —0.89 ±0.44 
and = —1.62 ± 0.56. The distributions of the second spectral index also 
show different mean values for the two samples, i.e. = —0.17 ± 0.37 and 



The neutrino break energies depend on the luminosity, redshift and the boost 
factor of the burst. The redshift distribution is displayed in Fig. 4. The variabil- 
ity sample has a significant contribution in the distribution at higher redshifts 
(0.8 < \ogz < 1.5), compared to the lag sample. Both distributions have their 




0.63 ±0.62. 
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Fig. 2. Distribution of the first spectral 
index a^. Solid line: variability sample; 
Dashed line: lag sample. Mean values 
with la errors are indicated for each 
distribution. 



Fig. 3. Second spectral index, fiy, of 
the neutrino spectrum for bursts in the 
variability (solid line) and lag sample 
(dashed line). Mean values with la er- 
rors are indicated for each distribution. 
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Fig. 4. Redshift distribution of the sources, solid line: variability sample; dashed 
line: lag sample. 

maxima around log 2; ~ — 1 which is consistent with the recent cognition 
that GRBs seem to follow Star Formation Rate (SFR). The boost factor is set 
to a constant value of F = 300 as it has been discussed above. 



There resulting break energies are shown in Fig. 5 and 6. Both samples scatter 
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Fig. 5. Distribution of first break energy. 
Solid line: variability sample; Dashed 
line: lag sample. 
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Fig. 6. Distribution of second break 
energy. Solid line: variability sample; 
Dashed line: lag sample. 



widely in both break energies, so that the errors in the average values, = 

5.0 ± 0.8 and 4'"^ = 5.4 ± 0.5 as well as ef"^ = 7.0 ± 0.9 and e^^^^ = 7.3 ± 0.7 
allow a deviation of about up to an order of magnitude from the mean values. 
Scatter plots of e|j, and el also show that there is a relatively strong correlation 
between the first and the second break energy in both samples with oc el, 
see Fig. 7. 

Fig. 8 shows the normalization of the variability (solid line) and lag spectra 
(dashed line). The distribution of the variability sample shows a peak with a 
mean value at log[A'j,/GeV s~^ sr~^ cm~^] ~ —8.8 ± 0.6, while the maximum 
for the lag sample is slightly lower, log[AV/GeV s~^ sr~^ cm~^] ~ —9.0 ± 0.7. 
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Fig. 7. versus e^. A correlation be- 
tween the two break energies is present 
due to a correlation of Equ. 7 and 9. 
Black squares represent the correlation 
for the lag sample, red triangles repre- 
sent values for the variability sample. 
The black line shows et = ef,. 
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Fig. 8. Distribution of the normalization 
of the single burst neutrino spectra, 
solid line: variability sample; dashed 
line: lag sample. 
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Table 1 



Parameters for neutrinos flux calculations. The numbers quoted as "typical values" 
are to be taken as rough bench marks, since all of these parameters fluctuate strongly 
as emphasized in the text. The calculation of the luminosity distance at a redshift 
between z = 1 — 2 is done using cosmological parameters of $7^ = 0-3, = 0.7 
and h = 0.71. 
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variability 


-8.8 ±0.6 


-0.89 ±0.44 


-0.17 ±0.37 


5.0 ±0.8 


7.0 ±0.9 


lag 


-9.0 ±0.7 


-1.62 ±0.56 


0.63 ± 0.62 


5.4 ±0.5 


7.3 ±0.7 


amanda sub. 
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-0.93 ±0.45 


-0.24 ±0.37 


4.8 ±0.8 


6.8 ± 1.0 
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Table 2 



Mean neutrino spectra parameters for the three samples, variability (568 bursts), 
lag (292) and the variability subsample (82 bursts). The standard deviation to the 
mean values has been calculated as an error estimate. The values used by Waxman- 
Bahcall (WB) are given as reference values. 
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Fig. 9. Comparison of coincidence and 
average spectrum (variability). 



Fig. 10. Comparison of coincidence and 
average spectrum (lag). 



The bursts that are most interesting in the presented sample in order to seek 
for an optimal spectrum for neutrino detection with a large volume neutrino 
Cherenkov telescope are those bursts with low first break energies - that means 
that the spectrum decreases with a relatively hard spectral index up to high 
energies, Ejy > 10^. Furthermore, high normalization bursts, > will have 
a strong influence on the shape of the coincidence spectrum. The coincidence 
spectra are being discussed in section 3. 



2.2 Sub sample from variability for AMANDA data analysis 



The last four years of BATSE operation time, 1997-2000, coincides with the 
beginning years of the AMANDA experiment. That implies the possibility of 
an examination of the neutrino signal from a BATSE subsample. 105 BATSE 
bursts lie within AMANDA's field of view (~ northern hemisphere) as is 
examined in Since the lag sample only contains bursts from before 1997, 
it cannot be used to examine any AMANDA burst. The variability sample 
however contains 82 of the 105 bursts. The parameter distributions of these 
bursts will not be discussed in more detail, but they are displayed in appendix 
A. 
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Fig. 11. Comparison of coincidence and average spectrum (AMANDA- variability 
subsample) . 



3 Between coincidence and average predictions 



For each burst i in a sample, a prediction of the prompt neutrino flux from 
this source can be made as described above. To analyze GRBs with a large 
volume neutrino telescope like AMANDA, a sample of bursts is analyzed to 
increase the signal rate which makes it interesting to look at a coincident 
spectrum of burst samples. The total flux $ is presented in the following form 
in this calculation, 

$ = ^^^. (13) 
n 



Here, n is the total number of bursts in the sample. 

In this section, the coincidence spectra of each examined sample will be de- 
scribed in more detail. Fig. 9 shows the expected energy spectrum for the 
variability sample, the solid line representing the coincidence spectrum, the 
dotted line displaying the average spectrum calculated from the mean param- 
eters as given in table 2. There are large differences between the two spectra, 
especially concerning the spectral indices at low and high energies. The differ- 
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Fig. 12. 292 bursts (lag, dotted line), 
568 bursts from variability method 
(solid line) and a subsample of 82 
bursts from the variability sample 
which can be analyzed with AMANDA 
(dot-dashed line). The variability and 
lag sample coincidence spectra differ 
clearly from each other and also from 
the average, diffuse Waxman/Bahcall 
prediction. This is a strong indi- 
cation that it is necessary to treat 
bursts individually in an analysis in 
order to optimize the signal expectation. 



Fig. 13. Summary of the mean spectra 
of all samples. Major differences can be 
seen between the variability (solid line) 
and lag sample (dotted line) , mostly due 
to differences in the second spectral in- 
dex and normalization. All calculated 
mean spectra differ significantly from 
the average diffuse prediction made by 
Waxman and Bahcall (dashed line) . The 
AMANDA subsample is represented by 
the dot-dashed line. 



ences between coincidence and average spectrum in the lag sample are rather 
in the normalization than in the shape of the spectra, see Fig. 10. The reason 
for these deviations lies in the contribution of individual bursts parameters far 
from the average. These bursts are not considered in average calculations at 
all, but can be responsible for a significant change in the spectral shape and 
normalization. The average and coincidence flux for the AMANDA subsample 
is shown in Fig. 11. Both versions of the energy spectra follow the larger sam- 
ple with smaller deviations, because the distributions of the subsample show 
a similar behavior as the whole sample (see appendix). What can also be seen 
is that for all samples both average and coincidence spectrum differ from the 
spectrum that is considered as the standard diffuse spectrum. The differences 
between diffuse and lag spectrum are quite small, but the spectrum of the 
variability sample has significantly different features compared to the diffuse 
spectrum. 

To compare the different samples to each other, the coincidence fluxes of all 
samples are shown in Fig. 12. The solid line shows the variability sample 
spectrum and the AMANDA subsample is represented by the dot-dashed line. 
The lag sample is given by the dashed line. The diffuse spectrum from (j^ is 
indicated by the dashed line as a comparison. It can be seen at first glance 
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that the two total samples show a very different behavior. Here, the difference 
from the average diffuse spectrum is of special interest, since this shows that 
coincidence spectra do not seem to be treatable as average diffuse spectra as 
easily. It can be seen that apart from the global normalization of the spectrum, 
the slope of especially the second spectral index is relevant for the overall 
normalization of the spectrum. A similar effect can be seen when comparing 
the average spectra of the samples to each other and to the diffuse spectrum, 
see Fig. 13. Here, however, the average variability spectrum agrees relatively 
well with the diffuse spectrum while there is quite a large difference between 
lag and diffuse prediction. 



4 Estimate of Cherenkov detection rate 

An estimate of the detection rate R can be given by folding the expected 
neutrino flux at Earth as it has been calculated above with the probability of 
the detection of the neutrino, 

hadow 

{e,E,)<!>,dE,. (14) 

Here, Py^i{Ey, -Emin) is the probability that a neutrino interacts with a nucleus 
to produce either a muon - / = /i - or an electromagnetic cascade - / =cascade 
- which is detectable in a large volume neutrino detector. It can be written as 

P.^i = Na I dEi^ri{EuE^^,,) (15) 

where Na is Avogado's constant, r; is the range of the produced muon - I = fi- 
or cascade - / =cascade - within detection range and da/dEi is the differential 
charged current cross section for u interactions. -Emin is the energy threshold 
of the detector for an event detection. The cross section is determined using 
the parton distribution functions given by |lfl[ l. where the model of is 
used. 

Pshadow is the probability that the neutrino is absorbed by the Earth. It is 
given as 

Pshadow exp i-xie)/x) . (16) 

The neutrino absorption length X{6) is dependent on the angle of the incom- 
ing neutrino towards the nadir, 9. It is determined by the distance that the 
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neutrino travels through Earth and the Earth's density using ()21l ) and the 
atmosphere's density of the US standard atmosphere modeL For a description 
of X(9) and A see for example (jj^ . 



The rate of neutrino induced muons per burst is displayed in Fig. 14, the solid 
line representing bursts from the variability sample, the dashed line showing 
the bursts from the lag sample. It can be seen that there is no burst yielding 
the probability of a whole event in a square kilometer detection array. This 
indicates that a source stacking method is useful for an analysis to get a higher 
significance for an actual detection. Furthermore, the fits from the variability 
sample yields more neutrinos on average than the ones from the lag method. 
As it has already been stated in (0), it can be confirmed in this analysis that 
there are a few bursts with rates above average which will dominate the flux 
from the sample. Bursts below average give a small contribution to the total 
signal. This is demonstrated in Fig. 15: The y-axis shows the total muon 
neutrino rate per sample. 



Here, imin is the lower summation index and Ri are the corresponding single 
burst rates. It is successively shifted to higher indices, Vin = 1, 2, n with 
n as the number of sources in the sample. The sources are arrange such as the 
lower rates have lower indices. For instance, the lowest rate in a sample gets 
the index i = 1, Ri, etc. Thus, by increasing iram, less and less sources are 
considered, starting by removing the least luminous ones. Fig. 15 shows Rtot 
versus Ri . . It can be seen that bursts with very low rates do not contribute 

''mm ^ 

significantly. Approaching the mean value of the distributions, however, the 
total rate starts to decrease rapidly. This behavior is seen in both samples. 
If the Waxman-Bahcall (WB) flux is used as the input spectrum (blue lines), 
the flux is constant at first as well, but plunges down extremely, because 
the rate distribution does not scatter very much as is seen in Fig. 14. The 
variability sample for instance shows a drop-off about an order of magnitude 
before the WB-fiux would drop of using the same sources. At this point, the 
flux has already decreased by one event where all events are still captured in 
the WB-scenario. This one event comes from the sources below average in the 
variability sample - the remaining 12.7 events result from the upper part of 
the spectrum. This shows that most of the contribution actually comes from 
bursts above average. Compared to the input of real parameters, the WB- 
scenario would give about 4.5 events less. The lag sample in contrast would 
only yield 3.6 events compared to 4.7 events in the WB-scenario. Mean and 
total values for the different scenarios are given in table 3. 
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Fig. 14. Muon neutrino rates for the single bursts using a threshold energy of 
-Emin = 100 GeV. The rates are given per burst and per square kilometer. On 
average, the neutrino flux rates determined in the lag sample are lower than 
the ones for the variability sample (log(i?|f^/km~^) = —2.73 it 0.94 compared to 
log(i?^'^''/km^'^) = —2.05 lb 0.67). The blue lines represent the same burst samples 
(solid: variability; dashed: lag) with the Waxman-Bahcall spectrum as input flux. 
All parameters except the angle are fixed as described above. 





mean rate [km ^] 


log total z^^rate [km ^] 


variability 


-|^Q-2.05±0.67 


13.8 


WB [var. sample] 


]^Q-1.86±0.23 


9.12 


lag 


-j^Q-2.73±0.94 


3.6 


WB [lag sample] 


-|^Q-1.86±0.23 


4.7 


AMANDA subsample 


-|^Q-2.05±0.59 


1.54 



Table 3 

Mean neutrino spectra parameters for the three samples, variability (568 bursts, 
single burst parameters and WB), lag (292, single burst parameters and WB) and 
the variability subsample of AMANDA bursts (82 bursts). The standard deviation 
to the mean values has been calculated as an error estimate. 



18 




Fig. 15. Total muon rate of the complete sample (variability: solid black line, lag: 
dashed black line). The lower summation limit is plotted on the x-axis. Only a 
few sources below average (red lines) contribute to the rate. As a comparison, the 
samples have been regarded using the Waxman-Bahcall spectrum as an input flux 
(blue lines), solid: variability; dashed: lag. 
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5 Conclusions and Outlook 



It is known since the detailed analysis of BATSE energy spectra that the 
parameters of the Band fits are not universal, but can differ significantly with 
the individual bursts. This implies that also the neutrino spectra of single 
sources have very different appearances. 

The challenging task in the future will be to get a more accurate descrip- 
tion of the redshift estimators, since until today, the relations are still derived 
from a very small sample of bursts with measured redshifts. The first redshift 
measurements of Swift give hope to a significant enlargement of the sample of 
measured redshifts in the near future which makes an improvement of the red- 
shift estimators possible. Since it is known that GRBs follow Star Formation 
Rate (SFR), a strong redshift evolution is expected. Therefore, the redshifts 
are an essential parameter in neutrino flux calculations and the currently avail- 
able redshift estimators are the most accurate way of dealing with redshifts 
so far. Other parameters like the spectral indices can be derived directly from 
the burst fits and therefore only include the errors due to uncertainties in the 
measurement of the photon spectrum. 

In this paper, a prediction of the neutrino flux from individual GRBs could be 
made using the model of pronipt neutrino emission that has been developed by 
(P). We follow the ansatz of (j^) by using individual parameters for each burst 
instead of using average parameters as it has been done by (0). It could be 
shown that it is quite important to look at ORB neutrino spectra individually 
to get an accurate description of a signal from these sources. This is of special 
interest for large volume neutrino telescopes as AMANDA, since an accurate 
description of the potential signal is necessary in order to optimize an analysis. 
This work shows that the use of an average spectrum for analysis purposes 
implies the danger of a misinterpretation of the results. Particularly, it could 
be shown that even the average spectra of a sample of bursts is likely to differ 
significantly from the diffuse prediction made by (j^. 
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A Parameter distribution of the AMANDA subsample (82 bursts) 
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Fig. A.l. a^: first spectral index of the Fig. A. 2. second spectral index of the 
neutrino spectrum. neutrino spectrum. 
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Fig. A. 3. First break energy. 



Fig. A. 4. Second break energy. 
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Fig. A. 5. Normalization of the spectrum. 



Fig. A. 6. Redshift distribution. 
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Fig. A. 7. Rate distribution for the 82 bursts in the AMANDA subsample. 
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